Global Minimizers for Free Energies of Subcritical Aggregation 
Equations with Degenerate Diffusion 



Jacob Bedrossian 



Abstract 

We prove the existence of non-trivial global minimizers of a class of free energies related to aggregation 
equations with degenerate diffusion on R d . Such equations arise in mathematical biology as models for 
organism group dynamics which account for competition between the tendency to aggregate into groups 
and nonlinear diffusion to avoid over-crowding. The existence of non-zero optimal free energy stationary 
solutions representing coherent groups in R d is therefore of interest. The primary contribution is the 
investigation of a notion of criticality associated with the global minimizer problem. The notion arises 
from the scaling of diffusive and aggregative forces as mass spreads and is shown to dictate the existence, 
and sometimes non-existence, of global minimizers. 

1 Introduction 

We consider global minimizers of the non-convex free energy 

T{u)= ( $(u)dx-\[[ u(x)u{y)K{x-v)dxd V :=S(u)-]-W(u) 1 (1) 

with u e L\(R d ) := {u G L 1 ^) : u > 0} such that \\u\\i := M over R d with d > 2. We are interested in 
determining for which choices of /C, $ and M there exist global minimizers with mass M. We restrict to the 
case where $ and K, are non-negative. We refer to S(u) as the entropy and W(u) as the interaction energy. 
Free energies such as ([T]) arise in the study of aggregation equations and Patlak-Keller-Segel models with 
degenerate diffusion [2[23li[TTJ[3[T7ll3[Ti[20l[Ta[21]. A typical example is 

u t + V • (uWK. *u) = Au m (2) 

for m > 1 and u(t) e L 1 + {R d ). These equations are formally a gradient flow for (p} with = -^zr[U m 

under the Euclidean Wasserstein distance [THl H E] , however as (P) is not displacement convex in the sense 
of [IS], the established theory does not apply. Numerical simulations indicate that for certain choices of K, 
and m, there are compactly supported stationary solutions to ^ which also appear to be attractors [25] H3] • 
The purpose of this work is to provide a rigorous analysis of the existence of these stationary solutions. 
Specifically, we give sufficient conditions for when there exist non-trivial global minimizers to (|TJ) , which are 
formally stationary solutions to @ with optimal free energy. We only consider the case K, radially symmetric 
non-increasing, which in particular implies that the non-local interaction is purely attractive. The minimizer 
problem without diffusion and more general /C has been considered elsewhere and is generally of a different 
flavor (see e.g. [H|9]). 

A typical approach to proving the existence of a global minimizer is to take a suitably strong limit of 
infimizing sequences. On bounded domains this may be carried through without issue, however on M d one 
must deal with the possibility that the sequence can lose some or all of the mass in the limit. To overcome 
this, we use concentration compactness arguments to establish tightness of infimizing sequences. Indeed, (p} 
was considered in the original work of Lions [14] . In [14] , the existence of minimizers was proved without the 
use of symmetry arguments in a number of cases using the concept of strict sub-additivity. Lions shows that 
this is effective for kernels with slow decay at infinity, but cases such as K, £ L 1 are not treated in detail. 
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Proposition [T] below contains the relevant results. Therefore, the primary contribution of this work is that 
it applies to kernels with fast decay. In contrast to [2], we rely on symmetrization arguments and a scaling 
analysis described below, as we could not verify strict sub-additivity in these cases. 

We consider K, non-negative, radially symmetric non-increasing such that K. G L^^°(R d ) n LP(R d \ i?i(0)) 
for some 1 < p < oo and 1 < p < oo. We define to* = (p + l)/p. Note that we are only considering the 
case m* > 1. The methods below allow us to define to* = 1 if K, G L°°, however, if K, has (for instance) a 
logarithmic singularity, the methods require us to define to* > 1, in contrast to [2J. Recall the following well 
known inequality which can be considered a generalization of the Hardy-Littlewood-Sobolev inequality: for 
all / G L 1 n L m * and S > 0, 

f(x)f{y)K{x - y)l Bs (\x - y\)dxdy < C \\K\ Bs ||^,« \\f\\\^ < C \\Kl Bi || L p.~ ||/||?- m * (3) 

2p-l 

where Cq is a constant that depends on p and d. Similar inequalities are used in [5] [21 [5J [T2]. Let J^f := 
{u G L\ n L m : ||u||i = M } and denote inf uG y M F{u) :— Im- Note that by considering a suitable sequence 
weak* converging to zero, necessarily Im < for all M > 0. We take <&(«,) to be strictly convex, non-negative 
and we assume /C, $ and M satisfy, 

lim inf 4^ > lc \\ICl Bs \\ LP ^M 2 - m \ (4) 

for some S > 0. Note that the quantity on the left hand side of course need not be finite. Condition (Q| 
implies that the problem is not supercritical in the sense of [5J and if it is critical, then the mass is assumed 
to be strictly less than the critical mass. The purpose of this assumption is to ensure — oo < Im and that 
sequences {u n } C with sup T(u n ) < oo have uniformly bounded entropy and L m norms (see for instance 
[14"] or [2]). However, we point out that kernels considered here are allowed to be more singular and have 
less decay than those considered in [2j. Moreover, we do not need regularity assumptions on JC. In [2], the 
primary purpose of the limitations was to ensure solutions to © were unique in L\_ n L°°. As in [Tl], we 
also require that lim z ^ ^{z)z~ 1 = 0. 

Before stating the main theorem, we describe the notion of criticality associated with ((T|) which dictates 
the main existence results in this work. The profile decomposition in [14] applied to sequences in 3^a/ with 
bounded L m norm suggests that the primary difficulties for proving Im is attained for some u* £ will be 
ensuring the mass of infinimizing sequences does not split apart or vanish. Naturally, this leads to considering 
the relative balance of the entropy and interaction energy as the mass of a sequence in 3^m spreads out. For 
simplicity, consider the case S(u) = m 1 _ 1 J u m dx with m > m* > 1, which arises from power-law degenerate 
diffusion, and suppose /C e L l . Let u(x) e L\ n L m and define u x (x) = X d u(Xx). Then, 

\dm-d \d p r i / <r — 7/ \ 

Hux) = ^—[MZ - y y J <*)<y) Yd K (-^\ dxdy. 

The key observation is that as A — > 0, the second term behaves like (A d /2)||/C|ji||w||2. If to > 2 then for 
sufficiently small A, J-{u\) < 0, whereas if m < 2, this is no longer true. The case to = 2 is in some sense 
critical, since the entropy and interaction energy scale the same as A — > 0, and the sign in the limit only 
depends on the value of ||/C||i. This scaling analysis is an important step (Lemma [lj to the proof of the main 
theorem below. We note that m > 2 are the exponents for which the equation |[2]) may be formally re-written 
as a regularized non-local interface problem [IT) . 

Theorem 1. Let m* > 1, K, non-trivial, M > and the hypotheses described above be satisfied. Suppose 
3 X, < x < oo such that 

<P(z)= X z 2 + o(z 2 ), z^0, (5) 
and additionally suppose that either of the following holds 
(i) x = 0. 

(ii) < x < oo and 2\ < ||/C||i < oo. 
Then Im < and there exists a radially symmetric non-increasing u* G [VW such that J~{u*) = Im- 
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and 



Lions in [T3] states results concerning the case when ((5|) does not hold for any x ( see (Hi) below). The 
contrast between the results shows how the scaling analysis manifests in the minimizer problem. As mentioned 
above, Lions also presents results which apply to kernels with slower decay at infinity (see (iv) below). Note 
that when kernels have slow decay, it is possible to prove the existence of minimizers in a wider variety of 
cases than Theorem Q] provides. The following statement is a consequence of Corollary II. 1 and Theorem II. 1 
in [H]: 

Proposition 1. (Lions \14V Suppose either of the following holds 

(Hi) The mass M is sufficiently large and there exists 1 < v < 2 such that for all t > 1 and z > 0, 
(iv) 3 a G (0, d) such that Vt > 1 and V£ G R d , 

Kffl > r a /c(0 (6) 

lim = 0. 

Then Im < and there exists a radially symmetric non-increasing u* G 3^m such that J~{u*) = Im- 

Remark 1. Due to the lack of convexity, to the author's knowledge, uniqueness is largely unresolved except 
when /C is the Newtonian potential or has similar special properties [101 113] . In critical cases, it is known to 
be not unique at the critical mass [5]. 

Remark 2. In applications, $ is often negative near zero [7]. However, for problems with degenerate diffusion, 
one can show S(u) > —A/, and the methods here apply simply by modifying (TTJ by an irrelevant constant 
depending on M. However, for problems with more general such as the Boltzmann entropy = it log u, 
the methods here would have to be modified. Similarly, modifications would have to be made to treat 
potentials JC which are unbounded from below, such as the logarithmic potential. 

We now point out that the condition \\K,\\i > 2\ in (ii) is close to sharp. 

Proposition 2. Let S(u) = and ||/C||i < 2. Then for all M > 0, Im = and there exists no non-zero 
global minimizers of the free energy. 

Proof. Recall from above that any global minimizer u* G will satisfy S(u*) = \\u*\\2 < oo. By Cauchy- 
Schwarz and Young's inequality for convolutions, 

*„2 If /• .. , ,. ... ., , {, IC\U 



F(u*) = \\u*\\i --J J u*(x)u*(y)IC(x - y)dxdy > [ 1 - 

This is clearly a contradiction unless = 0, as Im < for all M > 0. □ 

We now prove the main theorem. 

Proof. (Theorem [lj To prove (i) and (ii) we begin with the following lemma, which is a restatement of the 
scaling analysis discussed above. 

Lemma 1 (Scaling Lemma). Let (i) or (ii) hold. Then VM > 0, 3 4> G n with F{4>) < 0. 
Proof. Let <j) G fl 3^a/ and consider the mass-invariant scaling (j)\{x) — \ d <j){\x). Then for R > 0, 

F{<t>\) = S(cf)x) — / / (p(Xx)<p(Xy)IC(x - y)dxdy 

- d dx --JJ ^{y) Yd K, 

< X d (/ ^f^dx -\JJ <Kx)t>(y)y d >C (^) 1 Br (\x-y\) dxdy 
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By the (O and <f> g C^°, for all e > and A sufficiently small (depending on cf>) such that, 

\- d <5>{\U{\x))dx = J \- 2d $(\U(x))dx < xHWl +eU\\l 
Therefore, since /C <E Lj oc and 4> <E C^°, for all e we may pick A small such that, 



Therefore, if % = 0, we clearly have Im < 0. Moreover, if < x < °°j then the sign of the right hand side 
does not depend on A or <j), and is negative only if ||/C1b r ||i > 2%. By choosing R sufficiently large, this is 
equivalent to ||/C||i > 2\- □ 

Lemma [T] provides conditions under which Im < 0, but in general this is insufficient to imply the existence 
of non-trivial minimizers. Without strict sub-additivity, we cannot directly apply the results of [2] to attain 
assertions (i) and (ii). To recover, we will use a symmetrization argument and the following lemma, which 
in general is strictly weaker than sub-additivity. 

Lemma 2. Let (i) or (ii) hold and M\ > M2. Then Im x < Im 2 - 

Proof. Let u n £ 34/ 2 sucn tnat nm n-s-oo f(« n ) — > Im 2 - By Lemma [U there exists v 6 n 3^m such, that 
F{v) < and \\v\\i = Mi — M 2 . Without loss of generality, by Riesz symmetric decreasing rearrangement 
and K, radially symmetric non-increasing, we may take u n and v to be radially symmetric and non-increasing, 
since applying a symmetric rearrangement will only decrease the interaction energy and leaves the entropy 
unchanged [12]. Therefore Ve > 0, 3R n (e) > such that \\u n l R d\ BR ^ \\ p < e for 1 < p < 00 (note this does 
not imply any kind of tightness). Choose x n such that B-R n (x n ) D suppw = and let u n = u n {- — x n ). Define, 

z n (x) = v(x) + u n (x). (7) 

By v, u n > we have ||z n ||i = M\. Now by the approximately disjoint supports, 

F{z n ) = J ®(z n (x))dx ~\ J J z„(x)z n (y))C(x - y)dxdy 

< [ $(u n (x))dx + I $(«(x) + e)dx - ^ / j z n (x)z n (y)K.(x - y)dxdy. 

J J supp v J J 

Notice that since v £ C^°, by the mean value theorem, 

<&{v{x) + e)dx < / $(v)dx + £(e), 

1 supp v J 

such that liminf e _>o £(e) = 0. Therefore, by K. > 0, 

T{z n ) < T(u n ) + T(v) + £(e). 
Since F(v) < we may choose e sufhciently small to ensure lim inf n ^oo J-(z n ) < lim inf n^oc J-(u n ) = Im 2 ■ d 

We now prove the assertions (i) and (ii). Indeed, let u n G }>m be such that T(u n ) — > Im- As above, we 
may assume that u n is radially symmetric non- increasing. We now show {u n } has a convergent subsequence 
in the strong L 1 topology. 

We follow the approach of Theorem II. 1 in [14]. Following the work contained therein, tightness up to 
translation is established using the profile decomposition lemma (Lemma 1.1). Accordingly, there exists a 
subsequence of {u n }, not relabeled, such that one of the following three possibilities occurs, 

(i) Tight up to translation: 3 {y n } C K d for which Ve > 0, 3 i? > such that j M d\ BR r y ) u ndx < e. 

(ii) Vanishing: MR > 0, lim„->.oo sup yeRti $ BR r y s u n dx = 0. 
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(iii) Dichotomy: 3 { u n } , { u n} ; { v n} C L+, such that u n = u n + u n + v n with, for i S {1, 2}, u x n u 2 n = u l n v n = 
0, u l n ,v n < u n , lirrin^oo dist(suppu^,suppu2) = oo and lim^oc [[ti£||i = M - a, lim^oo ||u£||i = a 
and linin^oo ||u„||i = for some a, < a < M. 

Vanishing does not occur: 

Vanishing is ruled out by Im < 0. Indeed, Im < implies limn^oo W(u n ) > 0. Assume for contradiction 
that the subsequence (not relabeled) {u n } vanishes as n — > oo. Let q £ [m* /(2m* — 2),p), R > and by 
Holder's inequality, 



W(un) = / J u n (x)u n (y)JC(x - y)dxdy 

<\W n \\l q / {2q -i)\\ lcl B R - 1 \\ q + \ I u n (x)u n (y)K(x-y)dxdy + M 2 \\lCl v .d\ BR \\ 00 

J J Br 1 <\x-y\<R 

j_ 1 )||/C1b h _ 1 || 9 + \\K1 Br \b r -i\\oo / Un(x) / u n (y)dydx + M 2 \\)Cl Md \ BR \ 

J J\x-y\<R 



< IWllg/(2g- 



By interpolation, |wra||2<j/(2g-i) is uniformly bounded by 2q/(2q — 1) < to*, and since {u n } vanishes we may 
deduce, 

lim inf W(un) < ||/C1p , ||, + WKl^sJ^. 

As i? — > co, the last term vanishes since /C is radially symmetric non-increasing and the first vanishes by the 
dominated convergence theorem and /C S Lj . Therefore, we have deduced 

lim inf W(u n ) < 0, 

n— >oo 

which is a contradiction to /m < 0. 
Dichotomy does not occur: 

Although we do not have strict sub-additivity, we will take advantage of the weaker property, Lemma [21 along 
with radial symmetry, to rule out dichotomy. Suppose for contradiction that dichotomy occurs. By Riesz 
symmetric decreasing rearrangement, recall that u n is radially symmetric non-increasing. This together with 
the properties of the profile decomposition u n = u\ + u„ + v n implies one of u n or u n converges to zero in 
L°° . In particular, one of the sequences must vanish, which is the advantage of radial symmetry. Assume 
without loss of generality that u n —± 0. First note by the disjoint supports, 

S(u n )>S(u n ) + S(u 2 n ). (8) 

Then, 

W{u n ) = W(u n ) + W{ul) + W{v n ) 

+ / / u n u n^( x - y)dxdy + v n uljC(x - y)dxdy + v n u 2 n K(x - y)dxdy. 



By the © and interpolation, for any 5 > 0, 



w(u 2 n ) < + no 

J \B s ]]^ir«„i|oo 



\B e hhl\\\ 



< \\K1bAM\uI\C*- 1 + \\Kl*\ B ,h\\<\\Z>- ( 9 ) 



Similarly for any 6 > 0, 

< IIVnllJH/ClR^Jloo + ||iCl B ,|| 1 ^.oo||» n ||?-™*l|«„||™:, (10) 

and for i £ {1,2}, 

v n u l n JC{x~y)dxdy < ||u„||i ||<||i \\fCl U d\ Bs \\oc + \\ICl Bs \\ LP ,oo |K||}~ m /2 |KIC* /2 ||<||i _m /2 |Kllm* /2 - 

(11) 
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Finally, let 



d n = dist(supp it„, suppM„). Therefore, 



// 




(12) 



Putting the estimates (|5j)- (|12[) together with d n — > oo, JC radially symmetric non-increasing, along with 
linin-Hxj ||u^||oo = lim„^oo S(ul) = lim^oo \\v n \\i = and the uniform boundedness of ||ujj| m *, ||« n ||m* 
implies, 



This clearly contradicts Lemma [5] and rules out dichotomy, leaving only that there is a subsequence of {u n } 
which is tight up to translation. 
Conclusion of proof: 

Following [14] one may now prove without modification that tightness is sufficient to extract a subsequence 
(not relabeled) such that u n — > u* strongly in L 1 for some u* g y^i with .F(it*) = Im- Strong convergence 
is due to J-(u n ) — > Im and the strict convexity of $(u) (see [2]). This concludes the proof of assertions (i) 
and (ii). □ 
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